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Abstract 

Telomeres, ubiquitous and essential structures of eukaryotic chromosomes, are known to come in a variety of forms, but knowledge 
about their actual diversity and evolution across the whole phylogenetic breadth of the eukaryotic life remains fragmentary. To fill 
this gap, we employed a complex experimental approach to probe telomeric minisatellites in various phylogenetically diverse groups 
of algae. Our most remarkable results include the following findings: 1 ) algae of the streptophyte class Klebsormidiophyceae possess 
the Chlamydomonas-type telomeric repeat (TTTTAGGG) or, in at least one species, a novel TTTTAGG repeat, indicating an evolutionary 
transition from the /\raib/c/ops/s-type repeat (TTTAGGG) ancestral for Chloroplastida; 2) the /\raib/c/ops/s-type repeat is also present 
in telomeres of Xanthophyceae, in contrast to the presence of the human-type repeat (TTAGGG) in other ochrophytes studied, and 
of the photosynthetic alveolate Chromera velia, consistent with its phylogenetic position close to apicomplexans and dinoflagellates; 
3) glaucophytes and haptophytes exhibit the human-type repeat in their telomeres; and 4) ulvophytes and rhodophytes have unusual 
telomere structures recalcitrant to standard analysis. To obtain additional details on the distribution of different telomere types 
in eukaryotes, we performed in silico analyses of genomic data from major eukaryotic lineages, utilizing also genome assemblies from 
our on-going genome projects for representatives of three hitherto unsampled lineages (jakobids, malawimonads, and gon- 
iomonads). These analyses confirm the human-type repeat as the most common and possibly ancestral in eukaryotes, but alternative 
motifs replaced it along the phylogeny of diverse eukaryotic lineages, some of them several times independently. 
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Introduction 

The termini of linear eukaryotic chromosomes are protected 
by telomeres. Their DNA part is formed typically by a long 
array of conserved minisatellite sequences that tend to be 
conserved in particular groups of organisms, for example, TT 
AGGG in vertebrates and fungi (named here as the human- 
type; Meyne et al. 1989), TTTAGGG in most plants 
{Arabidopsis-type; Richards and Ausubel 1988), or TTAGG in 



insects (Okazaki et al. 1993; Frydrychova et al. 2004; Vitkova 
et al. 2005). These telomeres are maintained by a special re- 
verse transcriptase, telomerase, which elongates telomeres by 
addition of telomeric repeats and thus solves the so-called end 
replication problem (reviewed in Chan and Blackburn 2004). 
However, many exceptions to these rules are known, for 
example, telomeres maintained by retrotransposons in 
Drosophila melanogaster (Biessmann and Mason 2003), 
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diverse telonneric nninisatellite sequences in yeasts (Teixeira 
and Gilson 2005), plants with the hunnan-type or unknown 
telonneric sequences (Sykorova et al. 2003a, 2003b, 2006), or 
novel telonneric sequence in Arthropoda (Vitkova et al. 2005; 
Mravinac et al. 2011). An interesting diversity of telonneres 
was recently described also in the green algal group 
Chlannydonnonadales, where at least two evolutionary transi- 
tions fronn the ancestral TTTAGGG type to the TTTTAGGG 
{Chlamydomonas-type) occurred in the clade Chloronnonadi- 
nia and independently in a subclade of the Reinhardtinia 
clade; nnoreover, the hunnan-type telonneric sequence was 
found in sonne green algal species that fall within the Dunal- 
lielinia and Stephanosphaeria clades (Fulneckova et al. 2012). 

Eukaryotic algae are a polyphyletic assennblage of phyloge- 
netically diverse organisnns with different life styles and strat- 
egies and are thus of interest for telonnere biology, because 
they represent a substantial portion of the eukaryotic phylo- 
genetic diversity. Three algal lineages Glaucophyta, Rhodo- 
phyta, and Chloroplastida (the later including also land 
plants) represent direct descendants of an ancestral alga 
with a cyanobacteriunn-derived plastid and are thought to 
fornn a nnonophyletic "supergroup" called Archaeplastida, 
whereas other algal groups obtained their plastids fronn red 
or green algae through a process called secondary or tertiary 
endosynnbiosis (Archibald 2009). Different telonnere types 
were described in chronnosonnes of the nucleus and the 
nucleonnorph (vestigial endosynnbiont-derived nucleus) of 
cryptophytes or chlorarachniophytes, with the nucleonnorph 
telonneres presunnably descending fronn original telonneres of 
the ancestral-engulfed algal endosynnbiont (Gilson and Mc- 
Fadden 1995; Zauner et al. 2000). Previous genonne sequenc- 
ing projects reported hunnan-type telonneric sequences in the 
diatonns Thalassiosira pseudonana and Phaeodactylum tricor- 
nutum (Arnnbrust et al. 2004; Bowler et al. 2008) and an un- 
usual AATGe sequence in the red alga Cyanidioschyzon 
merolae (Nozaki et al. 2007). Synthesis of Arabidopsis-type 
telonneric repeats by telonnerase was observed in dinoflagel- 
lates (Fojtova et al. 2010; Zieike and Bodnar 2010), illustrating 
the diversity of telonneric sequences in Alveolata, where dil- 
ates possess telonneres with TTGGGG or TTTTGGGG repeats 
(Prescott 1994) and apiconnplexan taxa display several related 
telonnere types (see Kissinger and DeBarry 201 1 for review). 

Excavata is a potentially nnonophyletic "supergroup" of 
unicellular eukaryotes that ancestrally share a characteristic 
ventral feeding groove and an associated specifically orga- 
nized nnicrotubular cytoskeleton (Sinnpson 2003; Hannpl 
et al. 2009). Excavates nnay occupy a key position in the eu- 
karyotic phylogeny and include nnany innportant or biologically 
interesting species, yet they rennain annong the nnost poorly 
explored eukaryotic supergroups with regard to their nnolecu- 
lar genetic and genonnic features. This holds true also with 
regard to telonnere biology, as telonneres have been character- 
ized only in very few excavates, including the parasitic trypa- 
nosonnatids (Van der Ploeg et al. 1 984; Lira et al. 2007) and the 



diplonnonad Giardia intestinalis (=G. lamblia; Le Blancq et al. 
1991). Genonne sequences have been additionally reported 
for the parabasalid Trichomonas vaginalis (Carlton et al. 
2007) and the heterolobosean Naegleria gruberi (Fritz-Laylin 
et al. 201 0), but no infornnation is available on their telonneres. 
For other deep excavate lineages, such as Jakobida, 
Preaxostyla, or Malawinnonadida, no knowledge about telo- 
nneres and no representative genonne sequences are available. 

Here, we tested the activity of telonnerase and investigated 
the presence of nninisatellite repeats in diverse algal lineages 
using experinnental approaches. In addition, our on-going 
genonne sequencing projects enabled us to obtain the first 
data about telonneres in jakobids, nnalawinnonads, and gonio- 
nnonads (the latter group representing the closest heterotro- 
phic relatives of the cryptophyte algae). Finally, we surveyed 
diverse eukaryotic genonne sequences available in public 
databases and inferred their telonneric sequences by in silico 
analyses. This connbination of experinnental and bioinfornnatic 
analyses allowed us to describe the diversity of telonneres 
across the eukaryotic phylogeny (fig. 1), to confirnn a predonn- 
inant occurrence of the hunnan-type telonneric sequence in 
basal lineages, and to dennonstrate independent acquisi- 
tion of the sanne telonneric repeats in various phylogenetic 
lineages. 

Materials and Methods 

Algal Cultures and DNA Extraction 

The algal nnaterial used in this study originated fronn culture 
collections as specified in supplennentary table SI , Supplennen- 
tary Material online. Algae were grown in the reconnnnended 
liquid nnedia BBM or MASM (www.ccap.ac.uk/nnedia/pdfre 
cipes.htnn; last accessed February 21, 2013), or on nutrient 
agar plates. Clironnera velia was grown in a nnodified f/2 
nnediunn (www.ccap.ac.uk/nnedia/pdfrecipes.htnn; last access- 
ed February 21 , 201 3) in which natural seawater was replaced 
by seawater prepared by dissolving 23.38 g of "Red Sea" salt 
(Red Sea, USA) in 1 I of distilled water. Pliaeodactylunn tricor- 
nutum was cultivated in MASM nnediunn supplennented with 
30 nng of Na2Si03 • 9H2O per liter. 

The absence/presence of eukaryotic contanninants was 
nnonitored nnicroscopically and using algal cultures grown on 
BBM and bacterial LB agar plates. The identity of nnost of algal 
sannples was verified by polynnerase chain reaction (PGR) ann- 
plification (discussed later) and sequencing of the internal 
transcribed spacer and/or snnall subunit (SSU) ribosonnal 
DNA (rDNA) regions, and when necessary (in sonne green 
algae), phylogenetic analyses were conducted to confirnn 
the assignnnent into specific algal classes (data not shown). 
Genonnic DNA for PGR annplification was isolated using the 
"modified IRRI" method (Collard et al. 2007). Genomic DNA 
from a control alga Clilorella vulgaris (TEL01, supplementary 
table SI, Supplementary Material online) was isolated 
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Fig. 1. — Telomeres in a framework of the eukaryotic phylogeny. The schematic phylogenetic tree of eukaryotes has been drawn based on recent 
phylogenomic analyses (HampI et al. 2009; Parfrey et al. 2010; Burki et al. 2012; Derelle and Lang 2012; Laurin-Lemay et al. 2012; Zhao et al. 2012; Paps 
et al. 2013). Unresolved or contentious regions of the eukaryotic phylogeny are shown as polytomies; note that the supergroups "Excavata" and 
"Archaeplastida" are depicted as monophyletic for convenience, but monophyly of all their constituent lineages (specifically concerning the position of 
malawimonads and glaucophytes) remains uncertain. Lineages with no information available about their telomeres have been omitted for simplicity. 
Alternative, recently suggested positions of the root of the eukaryotic phylogeny are marked: root 1 (Rogozin et al. 2009); root 2 (Cavalier-Smith 2010); 
root 3 (Derelle and Lang 2012); root 4 (Katz et al. 2012). Telomere types documented for individual lineages are indicated on the left, with the various 
(putative) telomeric repeats indicated in different colors; NCT, noncanonical telomere (e.g., transposon based). Telomeric sequences confirmed by exper- 
iments or by a sufficiently complete genome sequence assembly are indicated in upper case, candidate telomeric sequences deduced from available draft 
genome sequences are indicated in lower case. When there is an apparently dominant telomere type known for a given taxon, the minor (typically 
secondarily derived) variants are shown in parentheses. For the sake of simplicity, the various unique telomeric sequences known in Saccharomycotina 
are indicated only by the collective label "various yeast repeats" and some of the additional unique fungal telomeric sequences are omitted (but presented in 
supplementary table S5, Supplementary Material online). Lineages that were targeted by experiments in this study yet failed to reveal their telomere structure 
are indicated with question marks. Details on the species representing the lineages in the phylogeny and the type of evidence for their telomere types are 
provided in table 1 and supplementary table S5, Supplementary Material online. 



according to protocol described by Saghai-Maroof et al. 
(1984). Isolation of DNA from other algal samples was per- 
formed according to a previously described protocol 
(Fulneckova et al. 2012), which involves the use of proteinase 
K during the lysis step. The DNA sample concentrations were 
estimated from agarose gels. 



Dot-Blot Hybridization 

Genomic DNA samples (~1 jig per sample) were dot-blotted 
onto Amersham Hybond-XL nylon membrane (GE Healthcare) 
and hybridized with radioactively end-labelled oligonucleotide 
probes (ATSB, CHSB, HUSB, TTCAGGG-SB, TTTAGGC-SB, 
T4AG2-SB, T3G3-SB, T2CG3-SB, Red alga-SB, supplementary 
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table S2, Supplementary Material online) as described in 
Sykorova et al. (2003b) with nninor nnodifications according 
to Neplechova et al. (2005). Briefly, membranes were hybrid- 
ized at 55 °C for 16 h and washed at 55 °C under low strin- 
gency conditions (2x saline sodium citrate [SSC] and 0.1 x 
sodium dodecylsulphate [SDS]); the final wash for the ATSB 
and HUSB oligonucleotides was done using a high-stringency 
washing buffer (0.6x SSC and 0.1 x SDS) to avoid 
cross-hybridization. Membranes for rehybridization with 
another probe were gently washed three times in 0.5% SDS 
at 80 °C. A control probe of mixed SSU and large subunit 
(LSU) rDNA fragments was prepared by mixing an equal 
amount of PGR products from several phylogenetically diverse 
algae (TEL213 Rhodella maculata, TEL97 Klebsormidium 
subtilissimum, TEL211 Tetraselmis chui, TEL207 Euglena 
geniculata, TEL01 Chi. vulgaris) obtained by amplification 
using a combination of gene-specific primers (18SrDNA-F 
and IBSrDNA-R [Katana et al. 2001] for TEL207, TEL211, 
TEL97, and TEL01; p4 and p23 [Van der Auwera et al. 
1994] for TEL97 and TEL211; and ITS-A [Blattner 1999] 
and p23 for TEL213); the mixture was labelled by DecaLabel 
DNA Labeling Kit (Fermentas, Thermo Scientific). The probe 
was used for final rehybridization overnight at 62 °G and 
low stringency conditions (2x SSG, 0.1 x SDS) or at 65°G 
and high stringency conditions (0.2 x SSG, 0.1% SDS). 
Membranes were exposed to autoradiography screens 
and signals were visualized using a phosphoimager 
FLA5000 (FujiFilm) and evaluated by the Multigauge software 
(FujiFilm). 

Restriction Digestion, Pulsed-Field Gel Electrophoresis, 
and Southern Hybridization 

Genomic DNA samples (1-5 |ig) were digested with restriction 
endonucleases Rsa\, Alu\, or Taq\ (NEB) and run on an 0.9% 
agarose gel in Tris-acetate-EDTA (TAE) buffer; the DNA frag- 
ments were alkali blotted and hybridized using the same 
hybridization and washing conditions as for dot-blots (for 
exceptions see supplementary material, Supplementary 
Material online). Agarose plugs with high-molecular-mass 
DNA samples for pulsed-field electrophoresis were prepared 
from lyophilized algal samples, and BAL-31 and restriction 
enzyme digestion was performed as described in Sykorova 
et al. (2006). Briefly, agarose plugs with high-molecular- 
weight DNA (TEL206 £ stellata, TEL207 £ geniculata, 
TEL133 Eustigmatos polyphem, TEL201 Vischeria punctata, 
TEL103 K. nitens, TEL97 K. subtilissimum, and TEL131 
Porphyridium purpureum) were digested with BAL-31 nucle- 
ase for 1 5 and 45 min (or 60 min), and then by the restriction 
endonucleases Sma\ (TEL97 and TEL131) or H/ndlll (TEL206, 
207, 133, 201, 103) (all enzymes from NEB). The DNA was 
then analyzed by pulsed-field gel electrophoresis using a GHEF 
Mapper (BioRad) under the following conditions: 1 % agarose 
(BioRad) gel in 0.5x TBE buffer, 6V/cm, pulses 0.5-35 s for 



20 h at 1 3 °G. Gels were alkali blotted and hybridized with the 
telomere probes. 

Telomere Repeat Amplification Protocol Assay 

Telomerase activity was investigated using a protocol originally 
developed for plant telomerases (Fitzgerald et al. 1996; 
Sykorova et al. 2003b) and applied with modifications to di- 
noflagellates (Fojtova et al. 201 0) and green algae (Fulneckova 
et al. 2012) (supplementary fig. S2A, Supplementary Material 
online). Briefly, 35-1 00 mg of lyophilized algal samples 
ground in liquid nitrogen were incubated in 1 ml of telomerase 
extraction buffer (Fitzgerald et al. 1996) and total proteins 
were recovered in supernatant ("crude extract") after centri- 
fugation at 17,000 x g for 15 min. The telomerase-enriched 
fraction was purified from the supernatant by precipitation 
with 10% polyethylene glycol (PEG) 8,000 and after centrifu- 
gation at 1 7,000 x g for 5 min the pellet was dissolved in one- 
quarter of the original volume of telomerase extraction buffer. 
Alternatively, the samples of crude protein extracts (without 
PEG precipitation) and the fraction of proteins not precipitated 
by PEG were used as specified in Results. A control telomerase 
extract was prepared from Arabidopsis thaliana Col-0 cell cul- 
ture. The amount of total protein in extracts was determined 
using the Bradford method (Bradford 1976). The telomere 
repeat amplification protocol (TRAP) assay was performed as 
described in Sykorova et al. (2003b) using a substrate primer 
47F (Fojtova et al. 2002) and a reverse primer TELPR30-3A 
(Fulneckova et al. 201 2) (supplementary material, Supplemen- 
tary Material online). Alternatively, different combinations of 
substrate primers, to cope with possible enzyme preference 
for the substrate primer sequence, and reverse primers repre- 
senting different telomere variants (supplementary table S2, 
Supplementary Material online) were used. Products were an- 
alyzed by polyacrylamide gel electrophoresis (PAGE), stained 
by GelStar(R) Nucleic Acid Gel Stain (LONZA) and visualized on 
a LAS3000 Imager (FujiFilm). TRAP products from selected 
algal species were cloned into the pCRIITOPO vector 
(Invitrogen) according to the manufacturer's recommenda- 
tions and sequenced (Macrogen). 

Gathering Genomic Data for Jakobids, Malawimonads, 
and Goniomonads 

Completed draft genome sequences and their systematic 
analyses of the jakobid Andalucia godoyi, the malawimonad 
Malawimonas califomiana, and the goniomonad Goniomonas 
avonlea will be published elsewhere together with details on 
DNA isolation, sequencing, and assembly protocols. Briefly, 
And. godoyi (ATCC PRA-185) and M. califomiana (ATCC 
50740) were sequenced by GS FLX Titanium platform (454 
Life Sciences/Roche) employing both shotgun and pair-end 
libraries. Draft assemblies were generated using Newbler 
2.6. Genome sequence data of the recently described species 
Gon. avonlea (Kim and Archibald 201 3) were generated using 
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the lllumina sequencing platform from multiple libraries, in- 
cluding two standard short insert libraries (300 bp) and mate 
pair libraries (2 and 6kbp) (Beijing Genomics Institute; McGill 
University and Genome Quebec Innovation Centre). Errors in 
the raw lllumina data were corrected using ALLPATHS-LG 
(Gnerre et al. 201 1), and the corrected reads were assembled 
using the ABySS de novo assembler (Simpson et al. 2009). 

Bioinformatic Analyses of Telomeric Sequences 

Candidate telomeric sequences were searched in our genome 
assemblies (discussed earlier) and in sequenced genomes 
available in various databases, including GenBank (www. 
ncbi.nih.gov; last accessed February 21, 2013), TraceArchive 
and Sequence Read Archive (SRA; http://www.ncbi.nlm.nih. 
gov/Traces/home/; last accessed February 21, 2013), Joint 
Genome Institute (www.jgi.doe.gov/; last accessed February 
21, 2013), Broad Institute (http://www.broad.mit.edu/; last 
accessed February 21, 2013), University of Tokyo (http://mer 
olae.biol.s. u-tokyo.ac.jp/; last accessed February 21, 2013), 
and EMBL (http://ct.bork.embl.de; last accessed February 21, 
201 3). Data sets of whole-genome sequence assemblies were 
downloaded in the FASTA format and searched using the 
BioEdit software and a string search for TTAGGG, TTTAGGG, 
TTGGGG, TTTTGGGG, TTTTAGGG, and TTAGG types of a telo- 
meric sequence. If unsuccessful, a search was then performed 
manually (by eye) at 5^- and 3^-ends of scaffolds looking for a 
TG/CA-rich telomere-like repetitive minisatellites. The position 
and distribution of candidate telomeric sequences (terminal 
and/or internal) in the genome assembly was subsequently 
assessed by an automatic search for the respective strings. In 
addition, unassembled genomic reads from Phaeocystis ant- 
arctica and Porphyra umbilicalis (Sanger reads in the TraceArc- 
hive and lllumina reads in the SRA archive, respectively) were 
searched using as a query a trimer of candidate sequences or 
one repeat of published telomere sequence of Cya. merolae 
(Nozaki etal. 2007). 

Results 

Sample Collection and Analyses of Algal Telomeres 

To cover a phylogenetically wide sample of algae, we culti- 
vated 48 algal strains from culture collections (supplementary 
table SI, Supplementary Material online), including members 
of Chloroplastida (Chlorophyta and Streptophyta), Rhodo- 
phyta, Glaucophyta, Haptophyta, Alveolata, Ochrophyta 
(Bacillariophyceae, Xanthophyceae, and Eustigmatophyceae), 
and Euglenozoa (supplementary table SI, Supplementary 
Material online). We examined all algal strains for telomerase 
activity by the TRAP assay (supplementary fig. S2/\, Supple- 
mentary Material online) and cloned the TRAP products from 
3 1 strains to determine what DNA sequence forms the ends of 
chromosomes (i.e., what sequence is synthesized by telome- 
rase) (table 1 and figs. 2 and 3; supplementary fig. S2, 



Supplementary Material online). Algal strains used in this 
study came both from groups where the telomeric sequence 
could be presumed from published data and from groups 
where the telomeric sequence has not been described yet. 
In the latter cases, a set of alternative reverse primer sequences 
were used in combination with three different substrate pri- 
mers to avoid false-negative results. In a subset of algal strains 
(32 in total), the occurrence of variant minisatellite telomeric 
repeats was examined by Southern hybridization (dot-blot hy- 
bridization and/or terminal restriction fragment [TRF] analysis; 
fig. 4 and supplementary fig. S3, Supplementary Material 
online) using telomeric oligonucleotide probes (supplementary 
table S2, Supplementary Material online). A terminal position 
of candidate telomeric sequences was tested by BAL31 nucle- 
ase digestion and Southern hybridization in eight algal species 
(representative samples shown in fig. 5, supplementary figs. 
S3 and S4, Supplementary Material online). 

Telomerase Activity Screening in Archaeplastida Using 
TRAP Assay 

We investigated 23 and 8 algal strains from Chlorophyta and 
Streptophyta, respectively, for the presence of a telomerase 
activity using the reverse primer TELPR30-3A with the 
Arabidopsis-type telomeric sequence. The Arabidopsis-type 
sequence was presumed as an ancestral telomere type for 
this group based on our previous results (Fulneckova et al. 
2012). Algal strains from the chlorophyte classes Chlorophy- 
ceae, Trebouxiophyceae, and Chlorodendrophyceae showed 
positive telomerase activity with products of a 7-nucleotide 
(nt) periodicity (fig. 2) and cloned TRAP products confirmed 
synthesis of the Arabidopsis-type telomeric sequence (table 1 ). 
A comparison of telomerase activity in the "crude" telome- 
rase extracts, the PEG-purified extracts, and the PEG-nonpre- 
cipitated protein fraction revealed that the telomerase activity 
is present also in the PEG-nonprecipitated fraction in all these 
telomerase-positive algal strains (supplementary table S3, Sup- 
plementary Material online). However, algal strains from the 
class Ulvophyceae failed to show a reproducible telomerase 
activity. Testing of reverse primers with sequences correspond- 
ing to alternative telomere types or minisatellite variants and/ 
or using different substrate primers to cope with possible tel- 
omerase substrate preference did not produce positive results 
(supplementary table S4, Supplementary Material online). A 
control experiment excluded the presence of telomerase in- 
hibitors in algal extracts (supplementary fig. SI, Supplemen- 
tary Material online; discussed later). In two cases, we 
experienced a very weak ladder of TRAP products using 
human-type reverse primer; however, we identified fungal 
contaminants in the respective two algal cultures by PCR 
(see Materials and Methods), which might be responsible for 
this residual activity in the samples tested. 

Three algal strains representing different branches of the 
streptophyte class Zygnematophyceae (TEL181 Zygnema 
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Fig. 2. — ^Telomerase activity in Archaeplastida investigated by TRAP assay. Telomerase activity in representative algal strains of Glaucophyta {A, TEL195 
G/a. nostochinearum), Rhodophyta {B, TEL213 R. maculata), Chlorophyta (C, TEL211 T. chui; D, TEL121 Dictyochloropsis irregularis; E, TEL94 
Pseudendodoniopsis botryoides and TEL124 Pseudendodonium basiliense), and Streptophyta (F, TEL198 IVIesotaenium endlidierianum; G, TEL97 
Klebsormidium subtilissimum, TEL100 K. dissectum, TEL101 K. flacddum, TEL103 K. nitens, TEL187 K. crenulatum) grouped according to their phylogenetic 
provenance (indicated above panels); the activity is shown using combinations of substrate and reverse primers — GG(21) and HUTC (human-type primer) {A), 
47F and TELPR30-3A (Arabidopsis-type) {B, C, F), or pSSyF and TELPR30-3A (D, E)- Synthesis of telomeric repeats corresponding to the human-type and the 
Arabidopsis-type sequence (compare with table 1) was observed in Glaucophyta (A) and three green algal classes (C, Chlorodendrophyceae; D, 
Trebouxiophyceae; F, Zygnematophyceae), respectively. Negative results were obtained in Rhodophyta (B) and Ulvophyceae (E). The Klebsormidiophyceae 
samples (G) showed synthesis of two different telomere types; alternative combinations of the substrate primer TS21 and the Clilamydomonas-type repeat 
reverse primer (CHTTRAPRevI), or of the substrate primer pSSyF and the TTTTAGG-type repeat reverse primer (T4AG2-C), displayed synthesis of a 7- or an 
8-nt periodicity of TRAP products (arrows) by telomerase of K. subtilissimum (TEL97) or other Klebsormidium spp., respectively (table 1). Differences in 
efficiency of telomerase purification during preparation from protein extracts are documented in samples shown on C, D, F, and G (summarized in 
supplementary table S3, Supplementary Material online). Triangles indicate different amounts of total protein (0.1 and 1 |ig) in protein extract without 
PEG precipitation (crude, cr.), in fractions nonprecipitated (supernatant, sup) and precipitated by PEG (telomerase extract, ex), except TEL94 {E: 0.1, 0.2 ^ig), 
TEL1 24 {E: 0.1,0.5 |ig), TEL97 (G: 0. 1 , 0.8 iig), TEL1 00 (G: 0. 1 , 0.3 iig), and TEL1 01 (G: 0. 1 , 0.4 iig). When one sample is indicated, 1 iig or a higher amount 
of total protein mentioned earlier was used, except TEL1 21 (D: all 0.5 |ig), TEL21 1 (C: sup 0.5 iig), TEL94 {E: cr. 0.1, sup 0.3 iig), and TEL101 (G: sup 0.4 iig). 
Telomerase-enriched extracts (50 ng of total protein) from Chlamydomonas hydra (TTTTAGGG), Arabidopsis thaliana seedlings (TTTAGGG), and Euglena 
stellata (TTAGGG) were used as a pattern control of an 8-, a 7-, and a 6-nt periodicity ladder, respectively; negative control (-), no extract. 
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Fig. 3. — ^Telomerase activity in algae outside Archaeplastida investigated by TRAP assay. Results of telomerase activity assay from representative samples 
of algal strains with telomerase synthesizing the human-type {A, Euglenophyceae, Haptophyta) and the Arabidopsis-t/pe {B, Xanthophyceae, Alveolata) 
telomeric sequence, and those with negative telomerase activity (C, Bacillariophyceae, Eustigmatophyceae), are shown. The ladder of positive TRAP products 
{A, B) corresponds to 6- or 7-nt periodicity of control samples (human- and Arabidopsis-type, respectively). The efficiency of telomerase purification during 
preparation in protein extract (summarized in supplementary table S3, Supplementary Material online) was monitored {A, B) without PEG precipitation 
(crude, cr.), and in fractions nonprecipitated (supernatant, sup) and precipitated by PEG (telomerase extract, ex), protein extracts containing lOOng and/or 
1 |ig of total protein were used. Negative results (0 were obtained using different amounts of total protein (indicated by triangle) from Phaeodactylum 
tricornutum (TEL231 ; 1 , 0.5, 0. 1 |ig), Vischeria punctata (TEL201 , 0.5 \xg), and Eustigmatos polyphem (TEL1 33, 1 , 0. 1 iig on left and 0.5 iig on middle panel) 
and different primer combinations (indicated under panels). Combinations of the substrate primer GG(21) and the human-type repeat reverse primer HUTC 
{A, Q or of the substrate primer 47F and the Arabidopsis-type repeat reverse primer TELPR30-3A {B, Q were used. Telomerase-enriched extracts (50 ng of 
total protein) from Chlannydonnonas hydra (TTTTAGGG), Arabidopsis thaliana seedlings (TTTAGGG), and Euglena stellata (TTAGGG) were used as a pattern 
control of an 8-, a 7-, and a 6-nt periodicity ladder, respectively; negative control (-), no extract. 



circumcarinatum, TEL196 Micrasterias crux-melitensis, and 
TEL1 98 Mesotaenium endlicherianum) displayed a positive tel- 
omerase activity and synthesis of the Arabidopsis-type telo- 
meric repeat (fig. 2 and table 1). In contrast, two different 
telomere types could be demonstrated in the class Klebsormi- 
diophyceae (fig. 2 and table 1). Four strains (TEL100 
K. dissectum, TEL101 K. flaccidum, TEL103 K. nitens, and 
TEL187 K. crenulatum) showed a pattern with 8-nt periodicity 
and cloning of TRAP products revealed synthesis of the Chla- 
mydomonas-type telomeric sequences. Interestingly, the TRAP 
assay failed with the Chlamydomonas-type reverse primer 
(fig. 2) in K. subtilissimum (TEL97) and showed a 7-nt period- 
icity that resulted from synthesis of a variant TTTTAGG minis- 
atellite repeat (table 1). Investigation of telomerase activity 
using various substrate primers showed a difference in sub- 
strate usage between telomerases from Zygnematophyceae, 
Chlorophyceae, Trebouxiophyceae, and Arabidopsis, resulting 
in a changed pattern of the TRAP products (supplementary 
fig. S2, Supplementary Material online). A similar difference in 
substrate primer usage was observed in plants from the 
monocotyledonous order Asparagales, which possesses telo- 
merase synthesizing human-type telomere repeats (Sykorova, 
Leitch, Fajkus 2006). 



Three algal strains covering different branches of rhodo- 
phytes failed to show telomerase activity when investigated 
using six variants of a reverse oligonucleotide primer derived 
from the telomeric sequence known in Cya. merolae (supple- 
mentary tables S2 and S4, Supplementary Material online) or 
four alternative reverse primers (supplementary table S4, 
Supplementary Material online) including the Arabidopsis 
type (fig. 2 and table 1). To check whether the failure of the 
TRAP assay could be caused by the presence of inhibitors in 
algal telomerase extracts, we performed a control experiment, 
in which the negative red algal telomerase extracts was added 
to a positive control extract from A. tlialiana in ratio 1:1 or 3: 1 . 
Neither of the red algal extracts posed a clear inhibitory effect 
(supplementary fig. SI, Supplementary Material online). 
Further analyses of telomerase activity using three different 
substrate primers with the human-type reverse primer 
showed a positive result in all three red algal samples, but 
only using the substrate primer 47F (supplementary fig. S5, 
Supplementary Material online). To verify the identity of the 
amplified minisatellite repeat, we performed TRAP reactions 
also with an alternative reverse primer (T3AG2-C) that was 
able to amplify successfully the TRAP products containing 
human-type repeats from other algal species (discussed 
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Fig. 4. — Dot-blot hybridization of genomic DNA with telomere and telomere-like minisatellite probes. Genomic DNA samples (1-2 |ig per dot) from algal 
strains and control samples (listed according to phylogeny position on left) were blotted and hybridized with radioactively labeled oligonucleotide probes 
representing different telomere types and derived sequences (indicated above dot columns). Samples in which the telomere type was revealed in a telomerase 
analysis (figs. 2 and 3, table 1) hybridized with the corresponding oligonucleotide, but occurrence of other minisatellites was also indicated. The T4AG2 and 
CHSB probes cross-hybridized and could not be distinguished by Southern hybridization (for details see supplementary material, Supplementary Material 
online). The T2CG3 minisatellite showed a signal across algal samples, but not in a terminal position (see supplementary material. Supplementary Material 
online); a similar situation was described in Aiiium (Sykorova et al. 2006). Control DNA samples represent the /Arab/c/ops/s-type {Chlorella vulgaris) and the 
human-type telomeres (human and Ipheion uniflorum, Alliaceae); note that the plant DNA contains also a portion of the ancestral Arabidopsis-type 
minisatellite. Control rehybridization of membranes was done with a mixed rDNA probe (see Materials and Methods for details); n.a., not analyzed. 
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Fig. 5.— Analysis of telomeres in Euglena stellata. Samples of 
high-molecular-weight DNA in agarose plugs were digested with BAL31 
nuclease for indicated times (in minutes) and then with restriction enzyme 
Hindi!!. DNA remaining in the plugs was analyzed by pulsed-field electro- 
phoresis (left panel) and the low molecular fraction of DNA after cleavage 
(diffused into reaction solution) was analyzed by conventional agarose gel 
electrophoresis (right panel). TRFs detected with a human-type telomere 
probe (without BAL31 digestion) range between 20 and 145kb and are 
sensitive to BAL31 digestion. A fraction of TRFs shortened after 15 and 
50min of digestion with BAL31 shifts to the low-molecular-weight frag- 
ments detected on conventional gel hybridization. The ladder is in 
kilobases. 

later) and could thus serve as a control. Sequencing of cloned 
TRAP products from these control reactions revealed that only 
repeats of the reverse primer sequence were amplified, sug- 
gesting a false-positive result. 

The glaucophyte alga Glaucocystis nostochinearum 
(TEL195) showed a telomerase activity (fig. 2) with a 6-nt pe- 
riodicity of the TRAP products using an Arabidopsis- or a 
human-type reverse primer. Cloning of the TRAP products 
from both primer combinations verified synthesis of the 
human-type telomeric sequence by the Gla. nostochinearum 
telomerase (table 1). Interestingly, algal strains from both 
streptophyte classes showed a substantial enrichment of tel- 
omerase in protein extract after PEG purification, in contrast to 
chlorophytes and the glaucophyte alga, which displayed sim- 
ilar activity both in purified and nonpurified samples (see sup- 
plementary table S3, Supplementary Material online). 

Telomerase Activity in Other Algal Groups 

We investigated telomerase activity also in representatives of 
other algal groups. According to published data, diatoms and 
euglenophytes should possess telomeres formed by the 
human type of telomeric minisatellite (Dooijes et al. 2000; 
Armbrust et al. 2004), but we were not able to detect any 



telomerase activity (fig. 3) in the diatom TEL231 P. tricornu- 
tum. In contrast, the three Euglena species tested (Eugleno- 
phyceae) and the haptophyte TEL2 1 0 Pavlova lutheri displayed 
high telomerase activity and synthesis of the human-type telo- 
meric repeats (table 1). Control experiments designed to in- 
vestigate possible preferences in substrate primer sequence or 
presence of inhibitors in diatom telomerase extract excluded 
these technical reasons of the TRAP assay failure (supplemen- 
tary table S4 and fig. SI , Supplementary Material online; refer 
the earlier discussion). An observed 7-nt periodicity and clon- 
ing of the TRAP products confirmed synthesis of the expected 
/\rab/c/ops/s-telomeric type in C. velia (fig. 3 and table 1), 
which is in agreement with its phylogenetic position within 
Alveolata close to dinoflagellates known from previous exper- 
iments to possess telomeres formed by the Arabidopsis-type 
sequence (Fojtova et al. 2010). Samples from the classes 
Xanthophyceae and Eustigmatophyceae showed very differ- 
ent results, despite the fact that they both belong to the algal 
phylum Ochrophyta within Stramenopiles. Although all five of 
xanthophytes (TEL95 Xanthonema cf. hormidioides, TEL202 
Pleurochloris meihngensis, TEL203 Xanthonema debile, 
TEL204 Heterococcus protonematoides, and TEL205 Botry- 
diopsis intercedens) showed telomerase synthesizing the Ara- 
bidopsis-type sequence, the two eustigmatophyte strains 
investigated (TEL133 Bus. polyphem and TEL201 V. punctata) 
did not reveal any reproducible telomerase activity (fig. 3). 
Similar to diatoms, control experiments using different com- 
binations of substrate and reverse primers showed negative 
result (supplementary table S4, Supplementary Material 
online), whereas a presence of inhibitors was excluded (dis- 
cussed earlier, supplementary fig. SI, Supplementary Material 
online). A PEG purification step was successful in telomerase 
enrichment of protein extracts from algal strains of Eugleno- 
phyceae, Haptophyta, Xanthophyceae, and C. velia (supple- 
mentary table S3, Supplementary Material online). 

Dot-Blot Hybridization Screening and Testing for 
Telomeric Localization of Minisatellite Repeats Using BAL 
31 Digestion 

We screened samples of algal genomic DNA by Southern 
hybridization using radioactively labeled oligonucleotides as 
probes (fig. 4) to unveil a possible occurrence of other 
telomere-like minisatellites in the respective genomes and to 
possibly identify candidate telomeric sequences in samples 
with no detected telomerase activity. We experienced difficul- 
ties in DNA extraction from several algal strains, mainly from 
Zygnematophyceae and rhodophytes, which showed the 
presence of colored substances and a low DNA yield; more- 
over, genomic DNA extraction was not successful for TEL196 
Micrasterias crux-melitensis and TEL198 Mesotaenium endli- 
cherianum. For the remaining samples from Chlorophyta, 
Streptophyta, Xanthophyceae, Euglenophyceae, Haptophyta, 
Glaucophyta, and C. velia, the dot-blot hybridization confirmed 
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the presence of teionnerjc minisatellites identified as "true" 
telonneric types synthesized by telomerase in the respective 
algal strains. However, dot-blot hybridization of genomic 
DNA fronn telomerase-negative strains did not suggest any 
other candidate telomeric sequence and in general, dot-blot 
hybridization signals were much weaker than we experienced 
in our previous study (Fulneckova et al. 201 2). A weak signal of 
control hybridization with a mixed probe consisting of a mix- 
ture of LSU and SSU rDNA sequences (see Material and 
Methods) may be caused by a wide phylogenetic span of our 
algal collection and a limited similarity among rDNA sequences 
or by a low quality of genomic DNA prepared by the proteinase 
K-based method, because we observed difficulties in PGR am- 
plification of control rDNA sequences and other Southern 
hybridization experiments (supplementary material. Supple- 
mentary Material online). The terminal position of a candidate 
human-type telomeric sequence in euglenophytes (TEL206 
and TEL207) and the terminal position of the Chlamydomonas 
type or the TTTTAGG type of a telomeric sequence in 
Klebsormidiophyceae (TEL103, TEL187, and TEL97) were veri- 
fied using BAL 31 nuclease digestion (fig. 5; supplementary 
figs. S3 and S4, Supplementary Material online). Subsequent 
rehybridization of BAL31 -digested samples of TEL97 {K. sub- 
tilissimum) with the Chlannydomonas-\)/pe sequence probe 
confirmed the presence of both sequence types in TRFs (sup- 
plementary fig. S3^, the bottom panel, Supplementary 
Material online). Investigation of the TRF lengths showed 
that the TTTTAGG-type sequences hybridize with 0.7-1. 5 kb 
long fragments (supplementary fig. S3A, right panel. 
Supplementary Material online), suggesting short telomeres 
similar to Chlamydomonadales. Correspondingly, digestion 
with Smal digestion produced longer restriction fragments 
and the signal of both TTTTAGG- and Chlamydomonas-type 
probes was distributed among multiple BAL31 -sensitive frag- 
ments of 2.5-23 kb length (supplementary fig. S3^, the 
bottom panel. Supplementary Material online). Besides these, 
short BAL31 -resistant fragments (1 .3-2.3 kb) representing in- 
terstitial telomeric sequences could also be seen in the hybrid- 
ization patterns of both probes. The presence of internal 
telomere repeats is also apparent in K. crenulatum (supplemen- 
tary fig. S4A, Supplementary Material online). Although high- 
molecular-weight restriction fragments hybridizing with 
Chlamydomonas-type probe shortened upon BAL31 treatment 
(supplementary fig. S4A the left panel, Supplementary Mate- 
rial online), the low-molecular-weight fragments were resistant 
to BAL31 , which reflects their internal (nontelomeric) positions 
(supplementary fig. S4/\, the right panel. Supplementary 
Material online). We also performed BAL31 digestion on 
both strains of Eustigmatophyceae to check whether the qual- 
ity of the genomic DNA could be the reason for the failure of 
dot-blot hybridization. Probing with the Arabidopsis-type or 
the human-type telomeric sequence, which are expected as 
candidate telomere types due to the phylogenetic position of 
Eustigmatophyceae in Ochrophyta, did not produce any 



specific signal (supplementary fig. S4, Supplementary 
Material online), thus confirming the negative results of the 
TRAP assay and dot-blot hybridization. BAL 31 nuclease diges- 
tion was performed also in TEL131 Porphyridium purpureum 
samples, but both investigated probes {Cyanidioschyzon-t/pe 
and human-type) failed to show any specific signal. 

Identification of Candidate Telomere Sequences in 
Genome Sequences of Phylogenetically Diverse 
Eukaryotes 

To cover a wider spectrum of phylogenetic lineages across 
the tree of eukaryotic life, we coupled our experimental in- 
vestigations with in silico searches for candidate telomeric 
sequences in published or publicly available genome se- 
quences, focusing on groups that have been ignored or 
poorly studied with regard to their telomeres. In addition, 
we take advantage of the genome data yielded by our 
on-going genome sequencing projects for three phylogenet- 
ically unique organisms, the jakobid And. godoyi, the mala- 
wimonad M. califomiana, and the goniomonad Gon. 
avonlea. We also used available genomic sequences to 
verify the presence of telomeric sequences that have been 
described previously for the respective organisms by meth- 
ods in telomere biology. We searched the genome assem- 
blies for stretches consisting of repeated units of the major 
known types of telomeric sequences (T^A^Go) and assessed 
them as candidate telomeric minisatellites by taking into ac- 
count their position and orientation with respect to adjacent 
sequences. Our simple database search could not uncover 
degenerated telomere types, like those known from yeasts, 
and experimental tests are also required to confirm the ter- 
minal position of the candidate sequences. We searched 143 
genomes (including 32 from species where the telomeric 
sequence has been published before) and 80 of them 
showed a convincing pattern of telomeric sequence (supple- 
mentary table S5, Supplementary Material online). A major- 
ity of the genomes that showed a dominant presence of the 
candidate sequence in terminal regions also exhibited inter- 
nal telomeric repeats occurring in short stretches or in large 
blocks (>100bp of uninterrupted minisatellite). The ge- 
nomes where we found only short or occasional repeats 
positioned terminally and/or in large internal blocks were 
considered inconclusive and ignored for the summary of 
the phyletic distribution of telomeric sequences in eukaryotes 
shown in figure 1 (except species with previously published 
telomeric sequences; supplementary table S5, Supple- 
mentary Material online). In several cases, we identified 
unexpected candidate telomeric repeats in genomes repre- 
senting hitherto unstudied key phylogenetic groups, for 
example, TTTCGGG in the parasitic relative of dinoflagellates 
Perkinsus marinus, TTTGGG in the heterolobosean 
N. gruberi, TTAGG in the labyrinthulid Aurantiochytrium 
limacinum, and the highly unusual 10-11 nucleotide 
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repeat unit TTTATT(T)AGGG in the rhodophyte Galdieria sul- 
phuraria (fig. 1 and supplementary table S5, Supplementary 
Material online). In addition, minisatellites that differed from 
the types "canonical" for the respective organismal groups 
were found in fungi and stramenopiles, indicating the evolu- 
tionary flexibility of the telomeric sequence at various phylo- 
genetic scales. Our database searches corroborated the 
experimental results from Haptophyta, Glaucophyta, and 
Chlorophyta, whereas no genome assemblies from Xantho- 
phyceae, Chromera, Ulvophyceae, Euglenophyceae, or dino- 
flagellates were available for analysis. The genome assemblies 
of two Nannochloropsis species (Eustigmatophyceae) dis- 
played the presence of large internal blocks of TTAGGG-type 
repeats in addition to several terminally positioned stretches; 
and without experimental evidence, these repeats should be 
taken as a candidate telomere sequence. Searches of unas- 
sembled genomic reads available for the red alga Por. umbili- 
calis did not identify the telomere types described in Cya. 
merolae genome or predicted in the Galdieria sulpliuraria 
genome (discussed earlier), but revealed a large number of 
reads containing a TTAGGG-type minisatellite. In several 
cases, these repeats could be assessed as internal, but whether 
any of the other sequences represent the true telomere cannot 
be verified without a full genome assembly or an experiment. 
A similar result was achieved when we searched unassembled 
Sanger reads from an on-going genome project for the hap- 
tophyte Plia. antarctica. 

Discussion 

People have long been fascinated by the question how mech- 
anisms of linear chromosome maintenance might have origi- 
nated. It is believed that recombination-based pathways, 
which today serve mostly as a backup mechanism (Fajkus 
et al. 2005), were original and were subsequently replaced 
by a more successful, steady, and efficient synthesis of telo- 
meres consisting of minisatellite sequences by telomerase. 
Telomeres of most investigated organismal groups are gener- 
ally conserved within the groups and conform to a limited 
number of minisatellite types. The repeat units of these minis- 
atellites are mainly variants of the T^A^Go sequence and their 
evolutionary success presumably depends on their properties 
stemming from the G-rich sequence and their capacity to form 
alternative DNA structures like the G-quartet or the T-loop 
typical for telomere function (de Lange 2005). The telomere 
sequence of the red alga Cya. merolae is somewhat atypical, 
but it is a G-rich sequence and genes for telomerase subunits 
were predicted in the genomic sequence (Nozaki et al. 2007). 
It is generally assumed that telomeres are formed by T/G-rich 
minisatellites and maintained by telomerase until alternative 
telomere and maintenance mechanisms are shown for a 
given organism. The results of our new experimental and in 
silico analyses have substantially expanded the sampling of 
telomere structures across the eukaryotic phylogeny. In 



addition, telomerase activity has also now been tested against 
a number of evolutionarily distant groups. The data enable us 
to paint a picture of telomere evolution across eukaryotes with 
an unprecedented level of detail (fig. 1), although many as- 
pects of the scheme remain to be clarified. 

Characterization of Telomeres in Many Algal Groups 
Remains Technically Challenging 

In algae, telomerase activity has been previously experimen- 
tally proven in photosynthetic and nonphotosynthetic dinofla- 
gellates (Fojtova et al. 2010; Zieike and Bodnar 2010) and in 
chlorophytes (Fulneckova et al. 2012). Despite a published 
completed genome sequence including predicted telomere 
sequences and a telomerase gene (Bowler et al. 2008), our 
TRAP assay was unsuccessful in the diatom P. tricornutum. 
This experimental failure could be for some unknown techni- 
cal reason or for the telomerase being active only in specific 
developmental stages (e.g., during auxospore formation) that 
were not presented in our samples at a detectable level. 
Another possible explanation comes from experiments in 
the silkworm Bombyx mori, where telomeres are maintained 
by telomerase with a very low processivity and thus difficult to 
be detected (Sasaki and Fujiwara 2000). Active telomerase 
synthesizing the Arabidopsis-type telomeric sequences was 
demonstrated by our experiments in C. velia (Alveolata), 
Xanthophyceae, chlorophytes, and streptophyte algae. In con- 
trast, we demonstrated the human-type telomeric sequence 
synthesized by telomerase in euglenophytes, haptophytes 
{Pav. lutheri), and glaucophytes (G/a. nostochinearum). 

Unexpectedly, we did not detect telomerase activity in the 
classes Eustigmatophyceae (Ochrophyta) and Ulvophyceae 
(Chlorophyta), and also Southern hybridization of genomic 
DNA with other minisatellite telomeric probes failed to resolve 
telomeric sequences, as was observed also in rhodophytes. A 
low signal intensity by Southern hybridization could be caused 
by impurities in the genomic DNA, which in our case had to be 
isolated by the proteinase K-based method (see supplemen- 
tary material, Supplementary Material online). Control TRAP 
experiments excluded the presence of telomerase inhibitors in 
protein extracts of all telomerase-negative algal strains 
(a diatom, Eustigmatophyceae, Ulvophyceae, and Rhodo- 
phyta). Our use of a different substrate primers should accom- 
modate variable telomerase preferences for substrate 
sequence. Despite these controls, it remains possible that 
there are technical reasons that prevent the telomere motifs 
from being resolved. Remarkably, we found a similar behavior 
of telomerases during PEG purification from protein extracts in 
various algal groups, suggesting their similar biochemical 
properties. In most cases, telomerase was enriched in a puri- 
fied fraction, with the exception of protein extracts from chlor- 
ophytes and a glaucophyte alga, which showed telomerase 
activity also in the nonprecipitated fraction (supplementary 
table S3, Supplementary Material online). However, PEG 
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purification led to removal of connpounds inhibiting TRAP 
assay from these extracts (Sevcikova et al. 2013). 

We did not detect telomerase activity in three red algal 
representatives investigated for synthesis of Cyanidioschy- 
zon-type telomeric sequences. This may be due to specific 
problems with PCR-based TRAP assay (a difficult G-rich tem- 
plate of a candidate sequence) or to occurrence of different 
telomeric sequences, because the investigated red algal spe- 
cies are only distantly related to Cya. merolae. The behavior of 
the three rhodophyte strains investigated here does not seem 
to be a peculiarity of one particular lineage, because they 
represent three different deeply diverged rhodophyte classes 
(fig. 1 and supplementary table SI, Supplementary Material 
online). Searches for telomeric motifs in databases suggest 
that some rhodophytes might have the human-type telomeric 
sequence, because a great portion of genomic reads from Por. 
umbilicalis (representing an additional class Bangiophyceae) 
contain this sequence type. However, many of them could 
be scored as internal sequences, and other telomere-like min- 
isatellites also occur in the genomic reads, making it difficult to 
distinguish the "true" telomeric type without an experiment. 
Moreover, the assembled genome sequence of Gal. sulphur- 
aria (a relative of Cya. merolae also belonging to 
Cyanidiophyceae) showed the presence of an unusual minis- 
atellite repeat (TTTATTAGGG or TTTATTTAGGG) predomi- 
nantly at the ends of scaffolds, suggesting a putative 
telomeric position. The sequence of this repeat unit seems 
to be an AT-rich variant derived from the typical T^A^Go telo- 
meric minisatellite. Neither the human-type nor other investi- 
gated telomere types (including the Cyanidioschyzon-t/pe) 
were found in the Galdieria genome assembly. The path of 
the telomere evolution in rhodophytes thus remains unclear. 
Assuming that Archaeplastida are monophyletic (but see, e.g., 
Burki et al. 2012) and that the human-type telomeric repeat 
occurs in the presumably basal archaeplastid lineage 
Glaucophyta as well as in most other major eukaryotic line- 
ages (discussed later), it is possible that the last common an- 
cestor of Archaeplastida had the TTAGGG telomeric 
sequence. This would then mean that this telomere type 
could be retained by some red algal lineages, but further test- 
ing of archaeplastid monophyly, and more evidence for the 
presence of the human-type telomere in rhodophytes, are 
needed to corroborate this scenario. 

Novel Telomere Forms Seem to Have Evolved in Some 
Algal Groups 

In plants and algae, two categories of evolutionary changes of 
telomere sequences have been described previously: 1) 
change of the telomeric sequence synthesized by telomerase 
to a related minisatellite sequence, that is, from the Arabidop- 
sis type to the human type in Asparagales (Sykorova et al. 
2003b) and from the Arabidopsis type to the Chlamydomonas 
type or the human type in Chlamydomonadales (Fulneckova 



et al. 2012); and 2) loss of a typical minisatellite and of telo- 
merase activity reported in the genus Allium (Asparagales) 
(Sykorova et al. 2006) and three genera of Solanaceae 
(Sykorova et al. 2003a). Our results bring evidence for further 
examples following the first and possibly also the second 
category. 

The first is exemplified by switches to variant minisatellite 
telomeric sequences in Klebsormidiophyceae, specifically a 
1-nt addition resulting in a change from the ancestral TTTA 
GGG {Arabidopsis type) to TTTTAGGG {Chlamydomonas type) 
early in the evolution of Klebsormidiophyceae, followed by a 
1-nt deletion in the lineage leading to TEL97 K. subtilissimum, 
which resulted in the unusual TTTTAGG sequence (this inter- 
pretation is based on the nested phylogenetic position of the 
latter K. subtilissimum within the genus Klebsormidium; Rindi 
et al. 201 1). The TEL97 strain still possesses a large amount of 
ancestral Chlamydomonas-t/pe sequence in its genome, 
which is similar to the situation in Hyacinthaceae (Asparagales) 
that have telomerase synthesizing the human-type sequence, 
but the ancestral telomere type still occurs in the genome at a 
level detectable by fluorescent in situ hybridization (Adams 
et al. 2001). The negative result of the telomerase assay 
using typical telomeric minisatellites as reverse primers in 
Eustigmatophyceae and Ulvophyceae suggests a change of 
the second type, but a more detailed study is needed to pin- 
point the evolutionary transition in the telomere sequence 
and/or maintenance mechanism. Similarities between land 
plants and algae are also seen in different substrate primer 
usage (see supplementary material. Supplementary Material 
online) that was demonstrated in Asparagales, and in that 
case its phyletic pattern did not show a simple correspondence 
to the phylogeny of this group (Sykorova, Leitch, Fajkus 2006). 
It seems that the substrate primer usage (supplementary fig. 
S2, Supplementary Material online) is linked to enzyme prop- 
erties of telomerases in different species, which generally 
show a varying accuracy in minisatellite repeat synthesis (Fitz- 
gerald et al. 2001 ; Sykorova et al. 2003b). 

Telomere Evolution in the Context of the Phylogeny of 
Eukaryotes 

Although the telomere structure remains unknown or enig- 
matic for many deeply branching eukaryotic lineages, includ- 
ing those with reference genome sequences available (e.g., 
Parabasalia, Archaemobae, or Apusomonadida), a combina- 
tion of previously published data with the results of our ex- 
periment and in silico searches allows to paint at least the 
major outlines of telomere evolution in eukaryotes (fig. 1). 
Metazoa (Traut et al. 2007; Gomes et al. 2010) and Fungi 
(Teixeira and Gilson 2005) show the TTAGGG-type as ances- 
tral to both groups and all changes in their telomeric se- 
quences seem to have happened later in their evolution. The 
same evolutionary path, that is, from the ancestral TTAGGG 
type to secondarily derived alternative telomeric repeats, can 



Genome Biol. Evol. 5(3):468^83. doi:10.1093/gbe/evt019 Advance Access publication February 9, 2013 



481 



Fulneckova et al. 



GBE 



now be deduced for Amoebozoa, Excavata, Stramenopiles, 
and probably also Archaeplastida (fig. 1). The TTAGGG 
motif remains the only telomeric repeat known in Choanofla- 
gellata, Haptophyta, and Rhizaria (excluding nucleomorph ge- 
nomes); hence, it is likely that it is also ancestral for these 
groups, but a much better sampling (especially for Rhizaria) 
is needed to confirm this. On the other hand, the plant type of 
telomeric repeats (TTTAGGG) may be ancestral in cryptomo- 
nads, which include cryptophyte algae and their heterotrophic 
sister lineage goniomonads. Inference on the type of the telo- 
meric repeat is difficult for the last common ancestor of Alveo- 
lata, because the two principal lineages, that is, Myzozoa 
(including dinoflagellates, perkinsids, apicomplexans, and 
chromerids) and ciliates, have different telomeres. 

Considering the wide occurrence of the TTAGGG telomeric 
repeat and its inferred ancestral presence in most major eu- 
karyotic groups, it is tempting to suggest that the human-type 
telomeric repeat was ancestral for eukaryotes as a whole. 
Deducing the ancestral state of any character for any taxon 
requires knowledge on the position of the taxon's root 
(the deepest branching point). The question about the root 
of the eukaryotic phylogeny has not yet settled, and at least 
four contradictory hypotheses, supported by different sources 
of evidence, have been suggested recently (fig. 1). Interest- 
ingly, all those root positions are compatible with the idea of 
the TTAGGG being the ancestral telomeric repeat for 
eukaryotes. 

Another important aspect to consider is the frequency of 
homoplasy in the evolution in telomeric sequences. Indeed, 
assuming the TTAGGG motif as ancestral, different phyloge- 
netic lineages evolved independently on each other to the 
same alternative motifs. A most notable case concerns the 
"plant" type TTTAGGG motif, which seems to have evolved 
independently in Chloroplastida, cryptomonads, fungi, oomy- 
cetes, xanthophytes, and alveolates (or at least their subgroup 
Myzozoa; fig. 1 and supplementary table S5, Supplementary 
Material online). Additional changes of secondarily evolved 
telomeric motifs then occurred many times, including shifts 
to novel motif (e.g., the TTTTAGGG motif in some chlamydo- 
monadalean and klebsormidiophyte green algae and in the 
apicomplexan Theileria [Sohanpal et al. 1 995], or the TTTCGG 
G motif in Per. marinus probably derived by a A-to-C substi- 
tution in the "plant" motif ancestral for Myzozoa) and rever- 
sions to the ancestral "human" type in some members of 
Chlamydomonadales (Fulneckova et al. 2012). 

In summary, our experiments and analyses substantially 
expand the sampling of telomere diversity in eukaryotes and 
strengthen the view of telomeres as evolutionarily flexible 
structures of eukaryotic genomes. The much more compre- 
hensive picture of the phyletic distribution of various telomere 
types should facilitate future studies on mechanistic causes of 
the evolutionary changes in telomeres and on underlying driv- 
ing forces. 



Supplementary Material 

Supplementary material, figures S1-S5, and tables S1-S5 are 
available at Genome Biology and Evolution online (http:// 
www.gbe.oxfordjournals.org/). 
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